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refraction  seismic  equipment  and  field  procedures  to  acquire  Rayleigh-wuve  data. 
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tions of  the  Rayleigh-wave  dispersion  technique  are  discussed. 
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FOREWORD 


This  report  is  « paper  that  was  presented  to  the  42nd  Annual 
International  Heating  of  the  Society  of  Exploration  Geophysicists,. 
The  Investigation  that  generated  the  Information  used  In  the  paper 
was  assigned  as  a work  unit,  "Techniques  for  Rapid  Subsurface 
Exploration,"  of  an  RDT6E  Project,  "Permanent  Construction  Materials 
and  Techniques  for  FY  70."  This  report  is  an  interim  report;  a 
final  report  will  be  published  in  FY  74. 

Engineers  of  the  U.  S.  Army  Engineer  Waterways  Experiment 
Station  (WES)  who  were  actively  engaged  In  the  investigation, 
analysis,  and  report  phases  of  this  study  were  Messrs.  R.  W.  Cunny, 
Z.  B.  Fry,  R.  F.  Ballard,  and  F.  K.  Chang,  of  the  Soils  and  Pave- 
ments laboratory  (S&PL) . The  work  was  performed  during  the  period 
December  1969  through  March  1972  under  the  general  supervision  of 
Mr.  J.  P.  Sale,  Chief,  StfL.  The  paper  was  prepared  by  Messrs. 
Ballard  and  Chang. 

00L  Levi  A,  Brown,  CE,  and  COL  Ernest  D,  Pelxotto,  CE,  were 
Directors  of  WES  during  the  conduct  of  the  investigation  and 
publication  of  this  report.  Mr.  F.  R.  Brown  was  Technical  Director. 
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CONVERSION  FACTORS,  BRITISH  TO  METRIC  UNITS  OF  MEASUREMENT 


British  units  of  measurement  used  in  this  report  can  be  converted 
to  metric  units  as  follows: 


Multiply 

feet 

pounds 

pounds  per  square  inch 
pounds  per  cubic  foot 
feet  per  second 


£* 

O.3OU0 
0. 145359237 
0.070307 

16.0185 

0.30U8 


Xa  Obtain 

meters 

kilograms 

kilograms  per  uquare 
centimeter 

kilograms  per  cubic 
centimeter 

meterB  per  second 
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Rapid  Subsurface  Exploration 

F.  K.  Chang*  and  R.  F.  Ballard,  Jr.**  j 

J 

This  paper  reports  an  investigation  of  improved  or  new  tech- 

i 

niques  to  extend  the  depth  and/or  resolution  capability  of  rapid 
shallow-depth  (less  than  500  ft)  explorations  that  would  quickly 

1 

and  economically  provide  an  engineer  with  accurate  information  on  1 

substrate  conditions.  Current  geophysical  techniques  for  rapidly  j 

exploring  construction  sites  and  investigating  foundations,  i.e., 
electrical  resistivity,  seismic  refraction,  and  vibratory  tech- 
niques, were  reviewed  for  relative  merits  and  shortcomings.  By 
combining  various  attributes  of  the  refraction  seismic  and 
vibratory  techniques,  it  was  considered  feasible  that  a study  of 
surface-wave  phenomena  could  measurably  extend  the  depth  limits 
of  current  investigation  techniques.  As  a result,  the  so-called 
"Rayleigh-wave  dispersion  technique"  was  envisioned  as  a potentially 
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promising,  in  situ  tent  method. 

A test  program  watt  formulated  to  adapt  conventional  refrac- 
tion seismic  equipment  and  field  procedures  to  acquire  Rayleigh- 
wave  data.  Theoretical  considerations,  test  procedures,  interpretation 
of  data,  and  comparisons  with  data  obtained  by  standard  methods  are 
presented.  Various  advantages  and  limitations  of  the  Raylcigh- 
wave  dispersion  technique  are  discussed. 
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INTRODUCTION 


Purpose  and  Scope 

The  purpose!  of  the  study  from  which  thiB  paper  was  compiled 
were  to  review  rapid  subsurface  exploration  techniques  currently 
employed  to  provide  engineering  information  on  subsurface  geo- 
logical formations  that  may  be  used  as  foundation  or  construction 
materials,  and  to  investigate  alternate  or  supplemental  techniques 
as  improved  tools  for  rapid  subsurface  exploration  (Ballard  and 
Chang,  1973).  The  scope  of  this  paper  involves  a cursory  review 
of  existini  rapid,  reliable,  low-cost  geophysical  methods  for 
site  investigation  and  location  oi  construction  materials,  as  well 
as  a detailed  description  of  a new  Rayleigh-wave  dispersion  tech- 
nique for  improvement  of  existing  rapid  exploration  capabilities. 

Background 

The  strength,  compressibility,  and  permeability  properties 
of  soils  or  rocks  usually  control  the  design  of  foundations. 

During  preliminary  site  or  foundation  investigations,  rapid  means 
of  indirectly  determining  these  properties  or  indexes  to  the 
properties  and  general  stratification  over  large  areas  are  desired. 
During  later  phases  of  the  investigations  at  particular  sites, 
more  detailed  or  specific  descriptions  of  the  foundation  materials 
are  desired.  The  subsurface  exploration  techniques  used  to  deter- 
mine the  desired  properties  of  the  material  car.  be  divided  into 
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Indirect  anvl  direct  techniques  (Glossop,  1968). 

Indirect  techniques 

In  geophysical  techniques  for  shallow-depth  exploration 
(maximum  depth  to  bedrock  » '.>00  ft),  the  seismic  refraction  and 
the  electrical  resistivity  techniques  are  commonly  used  (Dorbrin, 
1968,  and  Heiland,  19^0).  The  depth  as  well  as  the  velocity  of 
subsurface  materials  can  be  directly  determined  by  a seismic  refrac- 
tion sounding;  the  dynamic  confined  compression  modulus  of  elastic- 
ity can  be  calculated  from  the  velocity  if  the  density  of  the  mate- 
rial is  known  or  can  be  estimated.  The  electrical  resistivity 
technique  is  uced  to  detect  the  variation  in  resistivity  of  earth 
materials,  which  is  largely  dependent  upon  the  amount  and  salinity 
of  the  contained  water;  thus,  indirect  measiu  em  \itu  of  porosity, 
saturation,  and  permeubility  are  possible.  The  seismic  refrac- 
tion and  resistivity  techniques  each  have  certain  advantages  and 
disadvantages  as  exploration  tools.  Sometimes  the  two  techniques 
complement  each  other  for  certain  types  -f  subsurface  investi- 
gations, to  that  it  may  be  very  advantageous  to  use  both  as  a means 
of  achieving  speed,  economy,  and  reliability.  Ten  years  ago,  the 
U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES)  adopted 
a combination  of  a vibratory  technique  in  conjunction  with  the 
seismic  refraction  technique  to  determine  the  i\  situ  shear  and 
Young's  elastic  moduli  of  soils  (Fry,  1963).  The  vibratory  and 
seismic  refraction  techniques  were  used  to  measure  the  Rayleigh- 
and  compression-wave  velocities  of  soils,  respectively.  The 
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gravitational  and  magnetic  techniques,  generally  used  for  the 
study  of  regional  geological  structure,  are  without  depth  control 
and  have  little  significant  engineering  value. 

The  general  advantages  of  the  indirect  techniques  are: 

a.  Geophysical  techniques  of  exploration  identify  gross 
changes  in  character  of  subsurface  materials  by  sur- 
face measurement  of  the  changes  in  certain  physical 
properties  of  the  earth  at  or  near  its  surface. 

b.  Larger  areas  or  projects  of  greater  linear  extent  can  be 
explored  more  rapidly  and  economically  by  use  of  geo- 
physical techniques  than  by  use  of  borings. 

c.  The  surface  measurements  as  determined  by  use  of  these 
techniques  indicate  average  conditions  within  a limited 
area  and  not  along  a single  vertical  or  inclined  line, 
such  as  is  the  case  when  the  boring  technique  is  used. 

d.  Data  acquired  by  use  of  these  techniques  can  be  used 
to  detect  the  subsurface  irregularities  ,hat  often  are 
missed  by  borings. 

The  disadvantages  are: 

a.  Although  the  general  character  of  subsurface  materials 

can  often  be  estimated,  the  materials  cannot  be  definitely 
identified  by  geophysical  techniques  alone;  therefore, 
these  techniques  must  be  supplemented  by  borings  in 
which  representative  and/or  undisturbed  samples  are 
obtained. 
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b.  The  interpretation  of  geophysical  data  requires  trained 
and  experienced  personnel. 

Geophysical  techniques  are,  however,  especially  well  suited 
for  reconnaissance  and  preliminary  determination  of  formation  con- 
ditions of  large  areas  to  be  used  in  the  construction  of  missile 
sites,  military  underground  structures,  nuclear  reactors,  dams 
and  reservoirs,  tunnels,  highways,  airfields,  large  housing  pro- 
jects, etc.  Geophysical  techniques  have  also  been  used  success- 
fully for  the  determination  of  water  tables  and  deposits  of  gravel 
or  other  engineering  construction  materials. 

Direco  techniques 

The  conventional  or  direct  techniques  utilized  for  subsur- 
face exploration  and  determination  of  foundation  properties  include 
boring;  sampling,  measurements  of  penetration  resistance,  vane- 
shear  strength,  plate-bearing  capacity,  and  California  Bearing 
Ratio;  and  other  special  purpose  tests  and  the  conduct  of  desired 
laboratory  tests  on  disturbed  or  undj  sturbed  samples  as  the  situa- 
tion requires.  The  direct  techniques  of  exploration  are  relatively 
slow  and  expensive  and  are  not  discussed  further  in  this  paper. 

The  authors  feel  that  the  thrust  of  this  paper  should  be  directed 
toward  a description  of  the  Rayleigh-wave  dispersion  technique  as 
a potential  new  tool  for  rapid  subsurface  exploration  rather  than 
discussing  "yesterday's  newspaper.” 
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RAYLEIGH-WAVE  DISPERSION  TECHNIQUE 


It  is  generally  understood  that  the  conventional  seismic  re- 
fraction technique  will  provide  a compression-wave  velocity  and  an 
accurate  depth  determination  from  the  first  arrival  times  of  com- 
pression waves,  and  that  the  shear-wave  velocity  at  depth  can  be 
determined  by  the  vibratory  technique,  which  has  a limited  depth 
of  penetration  of  about  150-200  ft.  Shear  modulus,  Poisson's 
ratio,  and  Young's  modulus  can  be  calculated  by  using  the  compres- 
sion- and  shear-wave  velocities  and  the  mass  density  of  the  soil 
or  rock. 

Through  analyses  of  data  acquired  using  the  vibratory  technique, 
it  has  been  found  that  the  waves  produced  by  the  vibrator  are  pre- 
dominantly Rayleigh  waves.  It  is  considered  feasible  to  explore 
shallow  depth  by  means  of  Rayleigh  waves  generated  by  refraction 
seismic  tests  for  the  following  reasons: 

a.  Rayleigh  waves  as  well  as  compression  waves  are  generated 
by  the  refraction  seismic  tests. 

b.  Rayleigh  waves  are  more  easily  observed  than  shear 
waves. 

c.  A pulse  traveling  wich  Rayleigh-wave  velocity  should  be 
a prominent  feature  on  any  seismic  refraction  record 
unless  the  instrumentation  has  been  designed  specifically 
to  reject  it. 

d.  Rayleigh  waves  generated  by  earthquakes  have  been  used 
successfully  as  a tool  for  exploration  of  the  earth's 
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■j 


crust  and  mantle  structure  during  the  past  three  j 

decades.  ] 

In  view  of  the  reasons  mentioned  above,  a proposed  seismic  I 

a 

J 

refraction  test  procedure  and  a plan  of  test  to  record  the  com-  j 

1 

press  ion  wave,  shear  wave,  and  Rayleigh  wave  on  the  same  seis-  ,] 

I 

1 

mogram  have  been  considered.  Certain  preliminary  tests  have  been  i 

conducted  for  comparison  of  Rayleigh-wave  velocities  determined 

with  both  vibratory  and  refraction  seismic  techniques.  I 

Field  Procedure  1 

j 

The  field  procedure  for  recording  the  Rayleigh  or  R-waves  j 

during  refraction  seismic  tests  is  as  simple  as  that  for  record-  I 

ing  the  compression  or  P-waves.  Actually,  the  P-  and  R-waves  j 

can  be  recorded  on  the  same  seismogram  following  one  shot  if  the  < 

4 

instrumentation  incorporates  specially  designed  automatic  gain  \ 

J 

control  (AGC)  circuitry.  Considering  that  a sharply  breaking,  j 

large-amplitude,  first  arrival  P-wuve  is  necessary,  two  shots 
must  be  made  with  less  sophisticated  instrumentation:  one  using 

a large  amount  of  TNT  (3  to  5 lb)  for  recording  P-waves,  and  one 
using  a lesser  amount  of  TNT  (0.25  to  2 lb)  for  recording  R-waves. 

The  size  of  the  charge  depends  on  the  shot-to-detector  distance. 

At  a close  shot  distance  of  10  ft,  an  electrical  cap  may  be  suf-  , 

ficient.  I 

) 

For  routine  dynamic  foundation  studies,  a twelve-geophone  ] 

seismic  cable  with  25- ft  spacing  of  the  geophones  is  sufficient. 

i 

6 j 

I 

\ 


In  practice,  shots  are  fired  at  increasing  distances  from  each  end 
of  the  cable  spread  until  seismic  waves  from  the  bedrock  are  recorded 
on  an  oscillograph.  For  rapid  reconnaissance  in  a large  area,  a 
twelve-geophone  seismic  cable  with  geophones  at  50-ft  intervals 
can  be  used.  For  rapid  subsurface  exploration,  two  seismic  pro- 
files of  fairly  long  length  perpendicular  to  each  other  can  usually 
be  completed  within  two  days. 

The  instrumentation  system  employed  in  field  explorations 
is  similar  in  most  respects  to  that  used  in  shallow  refraction 
explorations.  The  band  of  frequency  response  including  the  geo- 
phone amplifier  and  galvanometer  is  wider  (around  the  range  of 
2 to  200  Hz)  than  in  standard  seismic  equipment.  The  geophones 
are  of  the  velocity  type  and  especially  designed  for  low-frequency 
measurement.  Present  WES  recording  seismic  instruments  were 
adapted  for  thv.  purpose  of  conducting  the  field  experiments. 

Data  Interpretation 


Definitions  of  the 
group  and  phase  velocities 

To  analyze  the  dispersion  characteristics  of  an  R-wave,  the 
group  or  phase  velocities  of  several  frequency  components  of  the 
R-wave  must  be  measured.  The  group  velocity  can  bo  measured  by 
assuming  that  the  wave  originates  as  a sharp  pulse  at  the  instant 
of  the  shot  and  at  the  surface  directly  above  the  shot  point.  On 
this  assumption.  Grant  and  West  (1965 ) determined  uhat  the  ^roup 


7 


velocity  vg  for  a travel  time  t of  a particular  cycle  of 
period  T observed  at  the  minimum  horizontal  distance  X from 
the  shot  is  given  by 


g 


X 

t 


(1) 


In  the  case  ox  soil  deposits,  the  seismic-wave  velocities 
increase  with  depth,  and  the  group  velocity  is  generally  less 

i 

than  the  phase  velocity.  The  phase  velocity  can  be  found  by 
numbering  the  peaks  and  troughs  of  the  oscillations  and  following 
them  from  one  adjacent  trace  to  the  nexc  as  will  later  be  de- 
scribed in  detail.  Assuming  that  the  period  of  a given  cycle  in 
the  sequence  does  not  change  drastically  from  one  trace  to  the 
next,  the  Rayleigh  phase  velocity  or  Rayleigh-wave  velocity  v 
can  be  approximated  by  the  relationship  given  by  Grant  an:’  West 
(1965), 

V - f , <2> 

r At 


where  AX  is  the  distance  between  adjacent  seismometers  and 
At  is  the  apparent  time  required  for  the  R-wave  to  travel  this 
distance.  (See  Figure  la  for  graphic  representation  of  the 
quantities. ) 

Seismogram  Analyses 

Figure  lb  is  a seismic  record  which  has  been  idealized  to 
show  the  prominent  wave  arrivals  at  two  adjacent  detector  loca- 
tions so  that  the  R-wave  analysis  could  be  better  illustrated. 
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GEOPHONES 

Ax 


DEFINITIONS  (IN  ORDER  OF  APPEARANCE  OR  USE): 

X = HORIZONTAL  DISTANCE  OF  SHOT  POINT  TO  FIRST  GEOPHONE  G1 

AX  = HORIZONTAL  DISTANCE  BETWEEN  GEOPHONES  G1  AND  G2 

♦ = RAYLEIGH-WAVE  ARRIVAL  TIME  (SUBSCRIPT  0 = FIRST  ARRIVAL, 
SUB-SUBSCRIPTS  I AND  2 = GEOPHONE  NUMBERS) 


APPARENT  TIME  REQUIRED  FOR  R-WAVE  TO  TRAVEL  AX  DISTANCE 
SUCCESSIVE  TROUGHS  OF  RAYLEIGH  WAVE  TRAIN 
AVERAGE  PERIOD 


At 

0,1 ,2,3,4 
T 

vr  : RAYLEIGH-WAVE  VELOCITY  OR  Vf  = 
L = WAVELENGTH  OR  L = VfT 

X 

Vg  = GROUP  VELOCITY  OR  vg  = ~ 


Ax 

At 


Fig.  1.  Layout  of  a tyj  .cal  seismic  line  and 
an  idealized  seismic  record 
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To  construct  R-wave  dispersion  curves  of  phase  and  group  velocity, 
the  following  steps  must  be  undertaken: 

Phase  velocity 

a.  Number  the  crests  or  troughs  of  the  Rayleigh-wave  train 
for  all  traces. 

b.  Read  time  of  arrival  for  each  crest  (or  trough)  and  the 
average  period  T for  each  trace. 

c.  Compute  the  incremental  time  At  between  adjacent  traces. 

d.  Determine  the  Rayleigh  phase  velocity  v for  each  wave 

r 

^ y 

number  by  = — where  AX  = distance  between  adjacent 
seismometers . 

e.  Compute  the  wavelength  L for  each  wave  number  by 

L ® v T. 
r 

f.  Tabulate  data. 

g.  Plot  v^  as  a function  of  L or  approximate  depth,  L/,?. 
Group  velocity 

a.  Number  the  crests  or  troughs  of  the  Rayleigh-wave  train 
for  at  least  one  trace. 

b.  Read  time  of  arrival  t for  each  wave  number  und  its 

o 

period. 

c.  Determine  group  velocity  v for  eaci  wave  number  by 

g 


ve  ” t 


(3) 


where  X * shot-to-detector  distance. 


Plot  v as  a function  of  T or  L (since  L = v T) 
g r 
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Data  presentation  and  reduction 


Providing  that  terrain  is  relatively  flat  and  background  noise 
levels  are  sufficiently  low,  the  oscillographic  data  presentation 
can  be  controlled  in  II.;  field  so  that  office  interpretation  time 
will  be  minimized.  This  can  be  accomplished  by  simply  adjusting 
the  galvanometer  alignment  so  that  each  trace  will  be  spaced  in 
proportion  to  the  distance  between  adjacent  geophones.  In  this 
ray,  the  interpreter  may  work  directly  from  the  field  records  if 
he  so  desires.  The  field  data  obtained  at  WES  test  site  2 (Figure  2) 
can  be  used  to  illustrate  this  and  several  other  interesting 
points.  As  will  be  observed,  data  were  recorded  in  four  225-ft, 
12-trace  se©ments  from  geophones  that  were  spaced  at  25-ft 
intervals.  The  size  of  the  seismic  impulse  source  ranged  from 
electrical  caps  to  1/3  lb  of  C-U  explosive  charges.  (The  size 
of  the  charge  was  kept  to  an  effective  minimum  so  as  not  to  deface 
the  landscape. ) As  a means  of  positive  R-wave  identification, 
traces  1,  2,  and  3 were  the  registrations  of  a transverse,  radial, 
and  vertical  geophone,  respectively,  located  at  the  first  station 
of  each  of  the  four  lines. 

Utilizing  the  simultaneous  relationship  of  the  vertical  and 
horizontal  (radial)  geophones,  one  can  plot  the  trajectory  of 
particle  motion.  Figure  3 shows  the  construction  of  these  dia- 
grams for  far-fieid  stations  located  from  250  to  700  ft  from  the 
seismic  source.  If  we  note  that  both  upward  vertical  motion  and 
outward  radial  motion  will  produce  a downward  deflection  of  the 
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•wave  dispersion  profile  (WES  test  site  2) 


UP  UP 


a.  250  FT  FROM  SOURCE 


UP 


b.  475  FT 
FROM  SOURCE 


UP 


c.  700  FT 
FROM  SOURCE 


Jp  LEGEND 

L 

DIRECTION  OF  PROPAGATION 

0.45  TIME  IN  SEC  DIAGRAM  WAS 
EXTRACTED  FROM  RECORD 

Fig.  3.  Particle-motion  diagrams  for  B-vaves 
(WES  test  site  2) 


seismogram  truce,  one  should  expect  oscillut ions  on  the  radial  trace 
to  lead  those  on  the  vertical  trace  by  90  degrees  in  phase.  ThiB 
is  seen  to  be  approximately  true  for  those  events  labeled  Rayleigh 
waves  on  the  record  shown  in  Figure  2.  Earlier  portions  of  the 
record  do  not  show  this  phase  shift  due  to  direct,  reflected,  and 
refracted  body  waves  causing  distortion  (Grant  and  West,  1965). 

In  each  of  the  three  cases,  the  particle  motion  roughly  takes  the 
shape  of  a retrograde  elliptical  orbit,  thus  fulfilling  the  defini- 
tion requirements  for  classical  R-wuves. 

Observing  the  first  or  near-fiold  record  shown  in  Figure  2, 
one  can  readily  see  that  the  trace  signature  characteristics  are 
completely  different  from  the  fur-field  records.  This  can  possibly 
be  explained  by  the  fact  that  the  distance  between  the  n..ot  and 
the  geophones  was  less  than  the  longest  wavelengths  that  were  ob- 
served at  the  far-field  stations;  therefore,  the  geophone  spread 
is  too  close  to  observe  dispersion.  However,  a second  factor  must 
also  be  considered;  the  seismic  source  for  the  near- li eld  record 
was  a single  electrical  cap  buried  only  O.S-ft  deep.  'Jhis  source 
very  likely  possessed  insufficient  energy  to  excite  the  lower 
frequencies  associated  with  depth  penetration  and  consequently 
greater  soil  masses. 

As  previously  stated,  interpretation  cun  be  performed 
directly  on  the  records  if  so  desired.  For  example,  consider  the 
second  record  in  Figure  2 (250-b75  ft)  for  illustrative  purposes. 
Beginning  with  trace  3,  some  11  distinct  oscillations  can  be 

lb 


observed,  numbered,  and  followed  through  the  succeeding  traces  on 
the  record.  A line  connecting  adjacent  troughs  can  then  be  con- 
structed in  a manner  analogous  to  a conventional  time-distance 
plot.  Essentially,  the  reciprocal  slope  of  this  line  will  yield 
the  wave  velocity.  An  average  period  should  then  be  determined  so 
that  the  wavelength  ca,»  be  computed  and  plotted  in  the  manner  pre- 
viously described. 

Conversion  of  wave  velocities  to  soil  properties 

To  become  a useful  tool  in  terms  of  conventional  soil  prop- 
erties, the  P-  and  shear-  or  S-wave  velocities  determined  by  the 
above  method  must  be  related  to  elastic  parameters  such  as  shear 
modulus  0 , Young's  modulus  E , and  Poisson's  ratio  v . This 
conversion  can  readily  be  accomplished  knowing  that  the  wave 
velocity  is  dependent  on  the  ratio  of  the  elasticity  of  the  medium 
to  its  mass  density  p and  on  the  wave  type.  The  relationship  of 

shear  modulus  G and  shear-wave  velocity  v and  density  p is 

8 

as  follows: 

G ■ v2  p (h) 

s 

S-wave  velocity  and  surface  F-vl ve  velocity  are  related  by  Poisson's 
ratio.  For  a homogeneous  medium  and  for  Poisson's  ratios  ranging 
from  0.2  to  0.5,  the  difference  in  velocities  is  less  than  9 percent. 
Therefore,  for  practical  purposes,  S-waves  can  be  considered  to  have 
the  same  velocity  as  R-vaves. 

Poisson's  ratio  v of  a soil  can  be  determined  by  the  rela- 
tion of  S-wave  velocities  v and  P-wave  velocities  v 

s c 
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Young's  modulus  E is  related  to  shear  modulus  G and  Poisson's 
ratio  of  the  soil  by 

E - 2(1  + v)  G (6) 

for  isotropic,  linearly  elastic  materials.  From  these  equations, 
values  of  E , G , and  v can  be  determined  through  the  measure- 
ment of  S-  and  P-vave  velocities,  provided  the  density  of  the  soil 
is  Known  or  can  be  estimated. 

Heukelom  and  Foster  (I960),  in  their  uynoinic  testing  of  pave- 
ments using  the  vibratory  technique,  found  that  excellent  correla- 
tions could  be  estimated  with  Known  depth  of  layers  if  they 
empirically  assumed  a measured  velocity  applicable  at.  a derth 
equal  to  half  its  associated  wavelength.  Uurir.g  the  past  ten  years, 
WES  has  extensively  employed  this  empirical  relationship  for 
investigations  where  boring  data  were  also  available  to  determine 
depths  of  interfaces  for  different  materials  and  agreement  has 
generally  been  good;  consequently,  this  relationship  is  assumed 
valid  for  the  R-wave  dispersion  method. 

Figure  U is  the  R-wave  velocity  for  WES  test  site  2 plotted 
as  a function  of  depth  (l/2)  prepared  using  the  R-wav^  dispersion 
technique.  In  addition,  data  obtained  by  using  the  surface  vibratory 
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technique  is  shown  for  comparative  purposes.  It  is  interesting  to 
note  that  excellent  agreement  exists  between  the  vibratory  and 
seismic  methods  to  depths  of  about  10  or  50  ft;  then  the  R-wave 
dispersion  data  diverges  to  indicate  higher  apparent  velocities. 

This  can  possibly  be  explained  by  the  fact  that  the  earliest  wave 
arrivals  (generally  th:  longest  wavelengths)  are  a mixture  of  re- 
flected F-,  S-,  and  R-waves  thut  tend  to  propagate  at  a higher 
apparent  velocity.  Unfortunately,  the  manual  method  of  determining 
the  phase  velocity  is  seldom  accurate  enough,  and  it  is  generally 
necessary  to  make  a Fourier  analysis  of  each  adjacent  pair  of 
traces  and  measure  At  from  the  phase  shift  of  the  Fourier  element 
whose  period  is  T (Grant  and  Wes*',  L'X>5  i • Cueh  on  analysis  can 
be  performed  with  available  computer  codes  if  it  is  considered 
necessary  by  the  data  user.  A second  possible  explanation  is  the 
difference  in  stress  levels  generated  by  the  t wo  methods  ; testing. 
Though  not  as  plausible  us  the  first  explanation,  il  must  be  con- 
sidered as  a possibility  and  further  study  must  bo  made  to  docu- 
ment possible  velocity  changes  as  a function  ci  ■ mrce  stress  level. 

Once  the  elastic  moduli  have  been  determined,  they  may  also 
be  conveniently  displayed  as  a function  of  dept1'..  !’hes<  lata,  which 
were  obtained  at  the  same  WES  test  site,  arc  sh'>'  n in  pigure  5. 

Once  these  parameters  are  supplied  to  a knowledgeable  designer,  a 
structure  with  predictable  performance  can  be  designed  (OCE,  1967)* 
Another  application  for  use  of  those  sell  parameters  voulri  be  for 
a dyne  .ic  analysis  of  existing  earthen  structures 
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such  as  dams 


potentially  exposed  to  earthquake  hazards. 

Theoretical  calculations 

As  a means  of  supplementing  the  interpretation  of  data  ob- 
tained by  the  R-wave  dispersion  technique,  comparisons  with  theo- 
retical dispersion  curves  based  on  soil  layering,  P-wave  velocity, 
S-wave  velocity,  and  density  of  the  media  can  be  made.  These  param- 
eters can  be  determined  by  borehole  sampling,  conventional  seis- 
mic refraction,  and  some  engineering  estimates  of  material  properties. 
Referring  again  to  Figure  U,  one  can  see  the  comparison  of  theoret- 
ical. and  experimental  values.  In  this  case,  all  three  methods  were 
used  as  an  aid  in  establishing  layer  thickness  and  properties  of 
the  Y/ES  test  site.  The  model  chosen  to  represent,  this  site  was 
representative  of  a four-layer  system.  The  fol  towing  characteristics 
were  assigned  to  each  succeeding  layer. 


Depth  Velocity,  fps  Wet  Density 


ft 

F-Wave 

S-Wave 

_ Pd/. 

0-12 

1,110 

350 

116 

13-56 

5,000 

650 

122 

57-76 

5,000 

9bo 

126 

77-“ 

10,000 

1*550 

Utilizing  these  assigned  descriptors  in  a R-wave  dispersion 
computer  code  developed  by  Watson  (1970),  the  theoretical  funda- 
mental mode  for  this  test,  which  was  treated  as  a four-layer 
problem,  can  be  calculated.  It  is  apparent  from  Figure  k that 
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some  differences  in  absolute  velocities  do  exist.  Future  research 


will  seek  to  reconcile  these  differences. 

^Evaluation 

The  advantages  and  limitations  of  the  R-wave  dispersion  tech- 
nique are: 

Advantages 

a.  R-waves  are  in  most  instances  easily  observed  and  should 
be  prominent  features  on  any  seismic  refraction  record, 

b.  The  computation  of  phase  velocity  of  an  R-wave  is  easy. 

c.  The  R-wave  dispersion  technique  is  a combination  of 
seismic  refraction  and  vibratory  techniques;  thus,  half 
of  the  field  working  time  could  be  saved  by  using  this 
technique.  However,  it  must  be  stated  that  additional 
interpretation  time  spent  in  the  office  tends  to  minimize 
the  overall  cost  savings. 

Limitations 

a.  The  investigation  area  must  be  large  enough  that  adequate 
information  can  be  obtained.  For  example,  if  the  bed- 
rock is  100  ft  deep,  the  length  of  the  profile  line 
(seismic-cable  spread)  should  not  be  less  than  500  ft. 

b.  Interpretation  of  the  R-wave  data  requires  the  services 
of  a competent,  trained  geophysicist. 

c.  The  R-wave  dispersion  technique  yields  velocities  that 
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are  consistently  higher  than  those  obtained  by  the  sur- 
face vibratory  technique  and  laboratory  test  methods. 


CONCLUSIONS 


Preliminary  field  experiments  indicate  that  the  R-wave 
dispersion  method  for  determining  shear-wave  velocity  by  surface 
measurements  is  indeed  feasible,  but  certain  questions  need  to  be 
resolved.  Most  apparent  is  the  lack  of  absolute  correlation  be- 
tween data  from  that  method  and  surface  vibratory  data.  At  cer- 
tain test  sites  examined  at  this  time,  the  R-wave  dispersion 
method  yielded  consistently  higher  velocities  at  depths  greater 
than  50  ft.  Indications  are  that  the  discrepancy  may  be  due  to 
wave  interference  (model  superposition  on  the  actual  seismograms) 
or  possibly  stress  conditions  at  the  source.  However,  additional 
study  is  needed  to  resolve  these  observed  differences. 
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